en 



o 

(N 






§ 
^ 



GEOMETRY OF QUANTUM DYNAMICS AND OPTIMAL 
CONTROL FOR MIXED STATES 

OLE ANDERSSON AND HOSHANG HEYDARI 



Abstract. Geometric effects make evolution time vary for different evolution 

curves that connect the same two quantum states. Thus, it is important to be 

able to control along which path a quantum state evolve to achieve maximal 

$H speed in quantum calculations. In this paper we establish fundamental rela- 

^1^ tions between Hamiltonian dynamics and Riemannian structures on the phase 

■^^ spaces of unitarily evolving finite-level quantum systems. In particular, we 

show that the Riemannian distance between two density operators equals the 
^_a^ infimum of the energy dispersions of all possible evolution curves connecting 

Cn the two density operators. This means, essentially, that the evolution time is 

a controllable quantity. The paper also contains two applied sections. First, 
we give a geometric derivation of the Mandelstam-Tamm estimate for the evo- 
lution time between two distinguishable mixed states. Secondly, we show how 
to equip the Hamiltonians acting on systems whose states are represented by 
invertible density operators with control parameters, and we formulate con- 
ditions for these that, when met, makes the Hamiltonians transport density 
operators along geodesies. 



^ 1. Introduction 

> 

Cn Geometric quantum mechanics is a branch of physics that has received much 

^^ attention lately. This is due in large part to the crucial role geometry plays in 

__ quantum information and quantum computing [TJ|ni[31lll|SJ|Sl[71[Hlini EH] • The 

performance of a quantum computer relies greatly on the efficiency of its algorithms, 

']^ and the ability to control which route an evolving state should take when joining 

2^ two given ones. This paper concerns fundamental aspects of the latter issue. 

,_^ A quantum system prepared in a pure state is usually modeled on a projective 

:• Hilbcrt space, and if the system is closed its state will evolve unitarily. Aharonov 

. ^^ and Anandan |llj showed that for unitary evolutions there is a geometric quantity 

^^ which, like Berry's phase [T21[T3], is independent of the particular Hamiltonian used 

H to transport a pure state along a given route. More precisely, they showed that 

the energy dispersion (i.e. 1/h times the path integral of the energy uncertainty) 

of an evolving state equals the Fubini-Study length of the curve traced out by the 

state. Using this, Aharonov and Anandan gave a new geometric interpretation of 

the time-energy uncertainty relation. 

Quantum computing and quantum information are theories developed for ma- 
nipulating mixed quantum states, i.e., statistical ensembles of pure states. Such 
states are usually represented by density operators. Many metrics on spaces of 
density operators have been invented to capture various physical, mathematical. 
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or information theoretical aspects of quantum mechanics. In this paper we make 
use of a construction inspired by Montgomery [14 to provide the spaces of isospec- 
tral density operators with Riemannian metrics, and we show that these metrics 
admit a generahzation of the energy dispersion result of Aharonov and Anandan 
to evolutions of quantum systems in mixed states. Specifically, we show that the 
energy dispersion of an evolving mixed state is bounded from below by the length 
of the curve traced out by the density operator of the state. Furthermore, we show 
that every curve of isospectral density operators is generated by a Hamiltonian for 
which the energy dispersion equals the curve's length. The latter result allows us to 
express the distance between two mixed states in terms of a measurable quantity, 
and we use it to derive a time-energy uncertainty principle for mixed states. 

Generically, mixed states of finite-level quantum systems are represented by in- 
vertible density operators. In a concluding section we define control parameters for 
Hamiltonians of systems whose states are represented by invertible density opera- 
tors, and give conditions for these that, when met, make the Hamiltonians transport 
density operators along geodesies. The section finishes with an application to al- 
most pure qubit systems. 

2. Geometry of unitary quantum dynamics 

In this paper we consider finite dimensional quantum systems in mixed states 
that evolve unitarily. They will be modeled on a Hilbert space T-L of dimension n, 
and their states will be represented by density operators. Evolving mixed states 
will be represented by curves of density operators, all of which, for convenience, are 
assumed to be defined on an unspecified domain < i < r. We write ViWj for the 
space of density operators on %. 

2.1. Riemannian structures on orbits of density operators. A density op- 
erator that evolves unitarily remains in a single orbit of the left conjugation action 
of the unitary group U{T-L) of T-L on ViT-L). The orbits of this action are in bijection 
with the set of possible spectra for density operators on T-L, where by the spec- 
trum of a density operator with fc-dimensional support we mean the nonincreasing 
sequence a = (pi,p2, • ■ • ,Pfc) of its, not necessarily distinct, positive eigenvalues. 
Throughout this paper we fix a, and write T^^a) for the corresponding orbit. 

To furnish I?(cr) with a geometry, let L(C'^,?^) be the space of linear maps 
from C'^ to T-L, equipped with the Hilbert-Schmidt inner product, and P{a) be the 
diagonal k y. k matrix with diagonal a. Set S{a) = {ip (L L{C'^,T-L) : V'^V' = P{<^)}j 
and define tt : S{a) — > T>{a) by 7r(V') = V'V'^- The fibration tt is a principal bundle 
with left acting gauge group U{cr) consisting of all unitaries in U{k) that commute 
with P{cr). (The action is C/ • -0 = ipU^ ■) Moreover, the real part of the Hilbert- 
Schmidt product restricts to a gauge invariant Riemannian metric on S{a): 

G{X, Y) = ^ Tr{X^Y + Y^X). 

We equip I?((t) with the unique metric g that makes n a Riemannian submersion. 
The Lie algebra if the gauge group is u((t). It consists of all antiHermitian k x k 
matrices that commute with P{(j). A connection for tt is given by the bundle of 
kernels of the mechanical connection form A on S{(j) defined by Aif, = II J,/,, 
where I^ : u(cr) -> u(cr)* and J^ : T^S{a) — > u((t)* are the locked inertia tensor 
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and moment map, respectively: 

I^^ 77 = G(VC^ W), MX) ■ i = G{X, i;S}). 

Vectors tangent to S{a) are called horizontal if they are annihilated by A, and a 
curve in S{a) is called horizontal if all of its velocity vectors are horizontal. Recall 
that for every curve p in 2?((t) and every ■i/'o in the fiber over p(0) there is a unique 
horizontal curve in S{a) that starts at ijjq and projects onto p. This curve is the 
horizontal lift of p extending from ■00. 

2.2. A geometric uncertainty estimate. Suppose A is an observable on T-L. Let 
Xa be the projection to 'D{a) of the gauge invariant vector field X^ on S{a) defined 
by 

We say that A is parallel at a density operator p if X^ is horizontal along the fiber 
over p, and we say that A is parallel along a curve p if A(t) is parallel at p(t) for 
every instant t. The precision to which the value of A can be known is quantified 
by its uncertainty function 



AA{p) = y Tr(yl2p) - Tr{Ap)^. 

In [IS] the current authors have proven that AA is bounded from below by h times 
the norm of the vector field Xa ■ 



(1) AA{p) > h^g{XA{p),XA{p)), 

(2) AA{p) = h^yg{XA{p),XA{p)) if A is parallel at p. 

The main argument is the following. For each ip in S{<t) there is a canonical 
identification between U{<j) and the fiber of tt containing tp, namely U i— >■ i/jW . The 
metric oiiU{a) obtained by pulling back G via this identification is independent of 
ijj- Restricted to u(cr) it is given by 

(3) i-v=lMi^^v + rik)P{<y))- 

Define the u(cr)-valued field £,a on ^{a) by tt*^^ — Ao X^, and write ^^ for ^a 
followed by projection onto the orthogonal complement of ilfe in u(cr). Then 

AA^ = h\g{XA,XA) + ii-ii). 

Now (IT]) follows from the observation that ^^ • Ca ^ 0, and (pi from the fact that 
^a{p) = if X^ is horizontal along the fiber over p. Note that for pure states, 
^^ = regardless of A since u(l) is spanned by ili. 

2.3. Distance, geodesies, and energy dispersion. The distance between two 
density operators with common spectrum a is defined as the infimum of the lengths 
of all curves in T>{(j) that connects them. There is at least one such curve whose 
length equals the distance, since V^a) is compact, and all such curves are geodesies. 
Moreover, horizontal lifting of curves is length preserving, tt being a Riemannian 
submersion, and a curve in V^a) is a geodesic if and only if its horizontal lift is 
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a geodesic in 5(cr), see [TB]. Here we show that the distance between two density 
operators pg and pi with common spectrum a satisfies 



(4) 



dist(po,Pi) = feinf / AH{p)dt, 
n- H Jo 



where the infimum is taken over all Hamiltonians H for which the following von 
Neumann equation is solvable: 

(5) p^hH,p] = XH{p), p{0) = po, p{t) = p,. 

in 

The length of a curve p in I?(o') is 

Length [p] = / ^Jg{p, p) dt. 

If p = Xh{p) for some Hamiltonian H, then, by (fl|, the length of /? is a lower 
bound for the energy dispersion: 



(6) 



Length[/9] < - AH{p)dt. 



There is a Hamiltonian H that generates a horizontal lift of p because the unitary 
group of H acts transitively on L(C'^,'H). For such a Hamiltonian we have equality 
in (|6| by (l2|. Moreover, we can take p to be a shortest geodesic. Then, 



dist(po,Pi)= \f AH{p)dt. 



Assertion Q follows. We refer to [TTj for a prescription how to produce a parallel 
transporting Hamiltonian from a given one, without affecting the evolution curve. 

2.4. A time-energy uncertainty relation. The Mandelstam-Tamm time-energy 
uncertainty relation [18 provide a limit on the speed of dynamical evolution. For 
systems prepared in pure states it implies that the minimum time it takes for a 
state to evolve to an orthogonal state is bounded from below by 7rS/2 times the 
inverse of the average energy uncertainty of the system. Recently, Jones and Kok 
[inilin] showed that the same inequality holds for mixed states, when orthogonality 
is replaced by distinguishability j^TJI^lj. Their proof involves an estimate of the 
rate of change of the statistical distance between density operators. Here we give 
a short geometric proof of this inequality. 

Consider a quantum system with Hamiltonian H, and suppose p is a solution to 
(pi). If po and pi are distinguishable, then 



(7) 



{AH)t>^, {AH) = - r AH{p)dt. 

^ T Jo 



To see this, let ipo ™ ^ "'^(Po) and ■01 hi 7r^^(pi) be such that dist(po,pi) = 
dist(i/)o, V'l)- The operators po and pi have orthogonal supports, being distinguish- 
able, and the same is true for ^o and ipi since the the support of ipo equals the 
support of Po, and likewise for -01 and pi. A compact way to express this is 

V-JVi = 0, vIV-o - 0. 
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If we consider ipQ and i/"! elements in the unit sphere in L(C , H), they are a distance 
of 7r/2 apart. In fact, V'(^) = cos{t)ipQ + sm{t)tpi, with domain < i < tt/2, is a 
length minimizing unit speed curve from ipo to ipi. Consequently, 

(8) dist(po,/Oi)>^/2. 

The relation (|7| now follows from ^ and (Is]) . Also note that the estimate (Isl cannot 
be improved. Direct computations yield ip'ip = P{a) and il^^tp = 0. Therefore, -0 is 
a horizontal curve in S{a), and hence ([8| is, in fact, an equality. From this it also 
follows that the estimate (|7| is saturated by any parallel Hamiltonian that generate 
p = ^-j/jt. 

3. Optimal Hamiltonians for mixed states of full rank 

Generically, the number of independent kets in a mixed state equals the dimen- 
sion of the Hilbert space. Such mixed states are represented by invertible density 
operators. From now on, we assume that the density operators in 'D(a) are invert- 
ible. 

3.1. Lie algebra controlled Hamiltonians. To achieve optimal computational 
speed in quantum computers it is desirable that the Hamiltonians transport states 
along shortest possible paths. Here we classify the Hamiltonians that transport an 
invertible density operator po along geodesies, and we provide conditions for control 
parameters of these Hamiltonians that, when satisfied, makes the Hamiltonians 
transport a ipo in the fiber over po along horizontal geodesies. 

As control space we choose the matrix Lie algebra u(n) , equipped with the metric 
given by ([3]). For each curve ^ in u{n) we define H^ by 

(9) H^ = ih^oP{a)-'/^ exp_ (J ^dt) (exp^ (J ^dt) P{<j)-^/^4. 

where P(cr)^^' ■^ is the diagonal nxn matrix whose j^ diagonal entry is 1/^/pJ, and 
exp_ is the negative time-ordered exponential. Also, let '0^ be the solution to the 
Schrodinger equation of H^ extending from ipQ. In the next section we show that 
every curve p extending from pq equals ip^Wfj for some curve ^ in u(n), and that 
ip^ is a geodesic if and only if ^ = adc ^, where ad J S, ■ f] — ^- [^,ri]. Furthermore, we 
show that ip^^ is horizontal if and only if ^ is contained in the orthogonal complement 
u(cr)-'- of u(o-) in u(n). 

3.2. Evolution operators that generate horizontal geodesies. The group 
U{'H) acts freely and transitively on S{(j) from the left, and the metric on U{'H) 
obtained by declaring the diffeomorphism U i-> UipQ an isometry is the left in- 
variant metric X ■ Y = Tv{{X^Y + Y^X)P{a)). An evolution curve UipQ is a 
geodesic in S{(j) if and only if [/ is a geodesic in IA{'H) extending from the identity 
operator. The second order geodesic equation in U{'H) can be reduced to a first 
order equation in 11(71) as follows. Define an isomorphism <j) : U{n) — >■ IA(T-L) by 
(j){U) = '4'oP{(T)^^/'^UP{a)~^/'^^\. Equip U{n) with the metric that makes cj) an 
isometry, and for each ^ in u{n) write Xi fo^' the vector field on UiTi) made up 
of left translates of d(p{£^). For a given curve U in l^iH) define a curve £, in u(n) 
by t/ = Xi{U). The Hamiltonian associated with ^, that transport ipo along C/f/'o 
in S{a), is H^ given by Q. Now, C/ is a geodesic if and only if ^ satisfies the 
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Arnold- Eulcr equation ^ = adt £,, see [531 and Appendix A Moreover, the inclusion 
of U{(j) in U{n) is an isometric embedding, when the former is equipped with the 
left invariant metric determined by ([3]). A straightforward verification shows that 
U^po is horizontal if and only if ^ is contained in u{a)-^, see Appendix [A| 

3.3. Geodesic orbit spaces and almost pure states. If a contains precisely 
two different, possibly degenerate, eigenvalues, every geodesic in ^{a) is generated 
by a time independent Hamiltonian. This since adt ^ = holds for every ^ in u(cr)^. 
To see this let rj be any element in u(ri), and write 



e = 





-4 



62 




v = 



mi 



■ni2 

'722. 



Then 



^■[Lri\^\ (piTr[^i24,r;n]+P2Tr[4Ci2,r]22]) =0 



since commutators have vanishing trace. (The corresponding result does not hold 
if CT contains at least three distinct eigenvalues.) Another way to put this is to say 
that 2?(<t) is a geodesic orbit space, i.e., a Riemannian homogeneous space in which 
each geodesic is an orbit of a one-parameter subgroup of its isometry group. 

By an almost pure state we mean a mixture of two pure quantum states in which 
one state is present in greater proportion than the other. Here we apply the above 
results to almost pure qubit systems. 

Two independent qubits are modeled by the standard basis elements ei and 62 
in C"^. Consider an ensemble of qubits prepared so that the proportion of qubits 
in state Cj is pj , where pi > P2- The initial state of the ensemble is represented 



by the density operator po — diag(pi,p2)- Chose ■00 = diag( 
over po, and let ^ be an arbitrary constant curve in u{a)'^: 

ee'- 

-ee-'^ 



'pi , Y^) in the fiber 



m 



e > 0, 0<t<l. 



The solution tp^ to the Schrodinger equation of H^ — ihS, that extends from ijjq is a 
horizontal geodesic, and the projection p = ^jV"! is a geodesic extending from po- 
Explicitly, 



P{t) 



Pi cos^ et + p2 sin^ et e*^ (p2 — Pi ) cos et sin et 

P2 cos^ et 



{p2 — Pi) cos st sin et pisin et- 

The curve p is a shortest geodesic between its end points provided that e is small 
enough, and dist(p(0), p(l)) = e. 



4. Conclusion 

The classic time-energy uncertainty relation by Mandelstam and Tamm implies 
that the evolution time between two distinguishable mixed states is bounded from 
below by a factor which is inversely proportional to the average energy uncertainty. 
In this paper we have shown that there is a fundamental relation between the 
length of evolution curves of mixed states, as measured by a specific Riemannian 
metric, and the energy dispersions of the Hamiltonians that generate the evolution 
curves. Our work thus indicates that the evolution time estimate derived from 
the Mandelstam- Tamm relation has a purely geometric origin. In fact, we have 
provided a geometric derivation of the same estimate. 
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In quantum computing it is desirable to have greatest possible control over the 
evolution of states. This to achieve maximum computational speed. Generically, 
mixed states of finite-level quantum systems are represented by invertible density 
operators. In the paper's concluding section, we have focused on quantum sys- 
tems whose states are represented by invertible density operators. There we have 
described how Hamiltonians acting on such systems can be equipped with control 
parameters, and we have provided conditions for these that, if met, guarantees that 
the Hamiltonians transport density operators along geodesies. 
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Appendix A. Constrained Arnold-Euler equation 

Suppose the density operators with spectrum a are invertible. Fix ipQ in S{<t), 
and identify U{n) and U{'H) according to the isomorphism (f) : U{n) — >■ U{Ti.) given 
by 

0(C/) = 4,oP{cj)'^'^UP{ct)-'/^4>1 

Then U{n) acts freely and transitively from the left on S{a). Thus U H> (l){U)ipo 
is a diffeomorphism from U{n) to S{a). Let Y^ be the push-forward of the one- 
parameter family of left invariant vector fields on U{n) generated by a curve ^ in 
u(n): 

We assert that the integral curves of Y^ are horizontal geodesies if and only if ^ is 
contained in u(cr)-'- and ^ = adt ^. The latter equation is called the Arnold-Euler 
equation [23]. 

The mechanical connection is such that a tangent vector is horizontal if and only 
if it is orthogonal to the fibers of tt. The tangent space at 0(C/)V'o of the fiber of it 
is spanned by the vectors (l){U)ipoV^ where 77 run through the matrices in u(cr). A 
straightforward computation yields 

G(re(0([/)^o),0(f/)W) = -^Tr((^t,^ + ^t^)p(^)). 

We conclude that Y^ is horizontal if and only if ^ is contained in u(ct)^. 
Let ip be an integral curve of Y^ . Its covariant derivative satisfies 

Moreover, by the Kozul formula [24, Prop 2.3], 

- XnG{X^,X^) - G{X^, [X^,Xn]) 
+ G{X^, [X,„X^]) + G{X„, [X^,X^]) 

= —2 adc ^ • 7] 

= — 2G'(Xad' {, Xrj) 
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for every 77 in u{n). Accordingly, S/^ip — X^^^, M). Thus, V is a geodesic if and 
only if ^ = ad J ^. 
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